Structure-dynamics relationship in coherent transport through disordered systems 
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Quantum transport is strongly influenced by interference with phase relations that depend sen- 
sitively on the scattering medium. Since even small changes in the geometry of the medium can 
turn constructive interference to destructive, a clear relation between structure and fast, efficient 
transport is difficult to identify. Here, we present a complex network analysis of quantum transport 
through disordered systems to elucidate the relationship between transport efficiency and structural 
organization. Evidence is provided for the emergence of structural classes with different geometries 
but similar high efficiency. In particular, we found a clear structural motif that renders transport 
properties robust against perturbations. Our results pave the way for a systematic rationalization 
of the design principles behind highly efficient transport which is of paramount importance for 
technological applications based on non-trivial quantum phenomena. 



Transport of charge or energy through disordered land- 
scapes is one of the most exciting mechanisms underly- 
ing biological and technological functionality . If the 
entities that are being transported behave wave-like, i.e. 
propagate coherently, interference resulting from scatter- 
ing off the disordered medium can result in strong focus- 
ing behavior due to constructive interference, as observed 
for example for an electron gas the coherent back- 
scattering of light from atomic clouds [5^ and predicted 
for Bose-Einstein condensates [6 . When this focus lies 
in the region to which an object should be transmitted, 
coherent behavior results in enhanced transport as com- 
pared to incoherent (particle like) processes [7, 8 . Con- 
sequently, it would be highly desirable to exploit such 
enhanced transport mechanisms. 

Constructive interference, however, relies on well- 
defined phase relations that need to be satisfied rather 
accurately, relations which get easily altered due to 
small changes in the geometry of the scattering medium. 
The onset of destructive interference then reduces trans- 
port, or suppresses it completely [9 . This results in a 
highly complicated structure-functionality relationship: 
two structures with hardly noticeable geometric differ- 
ences can lead to strongly different transport properties 
and two structures with similar transport potential might 
not share any geometric similarity. It is thus rather hard 
to identify geometric features associated with good trans- 
port properties, what would be absolutely necessary for 
the use of constructive interference as design principle in 
technological applications. 

The inherent complexity of the problem as well as the 
large number of degrees of freedom involved calls for the 
application of advanced statistical tools. Inspired by the 
substantial achievements of network science to elucidate 
complex systems like for instance economic growth [10], 
human diseases [11 and organic chemistry [12], our aim 
is to shed light on the elusive relationship between the 
structure of disordered media and constructive interfer- 
ence through the application of a network approach. 

Specifically, we considered a discrete two-level A^-body 



system, whose interactions are described by a tight- 
binding Hamiltonian 

N 3 

where cr^ describe the annihilation/creation of an ex- 
citation at site i, and J is the coupling constant. The 
interaction rate decays cubically with the distance be- 
tween the sites in accordance with dipole-dipole interac- 
tion. Within this model a structure is defined by the 
positions of the TV sites. The initially excited site (input) 
and the output site where the excitation has to be deliv- 
ered are located at diagonally opposite corners of a cube 
of side ro, while the remaining N — 2 sites are placed at 
random positions within this cube. 

While we looked at systems with A/" = 4, 6 and 8, 
most of our attention focused on A/" = 6, being the small- 
est set for which non-trivial behavior emerged. For this 
case, a large sample of 100 millions random structures 
was generated. This sample covered the whole spectrum 
of transport efficiencies where E is defined as the max- 
imal probability 

^ = max,e[o,r]|(m|e^^>iit)|' (2) 

to find the excitation at the output site \out) within a 
short time interval after initialization. In order to access 
fast transport that necessarily needs to result from con- 
structive interference, we chose r = -L ^'"'"3°'"* , i.e. a 
time-scale that is ten times shorter than the scale asso- 
ciated with direct interaction between input and output 
sites [T3] . 

Being interested in the characterization of efficient 
transport, our analysis focused on structures with E > 
0.9, a property that is satisfied by only 14280 configura- 
tions out of the generated 10^. From a structural point 
of view this reduced set is highly heterogeneous, meaning 
that two structures with similar efficiency do not nec- 
essarily share any evident common pattern [13 . Such 
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structural variety hinders a straightforward interpreta- 
tion of the geometrical features compatible with efficient 
transport. 

To unravel this structure-dynamics relationship, we 
applied a set of tools based on complex networks origi- 
nally developed for the characterization of molecular sys- 
tems im |TB| . Specifically, these methods are designed 
to analyze large ensembles of configurations, potentially 
allowing in this case a systematic classification of struc- 
tures which lead to exceptional transport. We gener- 
ated a complex network where the 14280 structures with 
^ > 0.9 represent the nodes and a link is placed be- 
tween them if two structures are geometrically similar 
independently on the specific dynamics of the excitation. 
The parameter used to estimate structural similarity is 
S'^ = YliLi ^i/^^ where di ^ i e {1, 6} is the distance 
between corresponding sites in any two structures. Links 
were placed when S'^ < 0.0125 Tq (see Methods in Supp. 
Mat. for further details on the network creation proto- 
col). The resulting quantum-structures network is shown 
in Fig. [l^. In this picture, nodes are proportionally close 
in space according to the amount of common neighbors. 
The color coding adopted for this network will be dis- 
cussed in detail below, anticipating at this stage that 
the presence of densely connected regions, i.e. clusters of 
nodes which are highly interconnected among each other, 
indicates the presence of groups of structures with com- 
mon geometrical motifs. 

We identified these regions using a network clusteri- 
zation algorithm based on a self-consistency criterion in 
terms of network random walks [1^ (see Methods in 
Supp. Mat.) that split the network into eight clusters 
comprised of structures with similar sites arrangements. 
The eight clusters have very different relative popula- 
tions, respectively of 73.3%, 9.2%, 4.7%, 3.2%, 2.6%, 
1.9%, 1.2%, 1.1%. As such, there is a bias towards cer- 
tain structure types with respect to others. Cumulatively 
the eight clusters represent the 97% of the whole sample 
with an unclassified 3% due to noise (see Fig. [s] in Supp. 
Mat.). 

The detection of eight clusters does not imply the pres- 
ence of eight distinct dynamical behaviors. We there- 
fore probed transport robustness against random dis- 
placements of the individual sites of a structure, with 
displacements restricted to a cube of side 0.05 ro cen- 
tered around the original position of the site. For each 
structure, AErand was calculated as the original effi- 
ciency minus the average efficiency obtained from 1000 
site-randomizations. In this scheme, structures were kept 
rigid under the assumption that the dynamics occurs on 
a much faster time scale than low-frequency fiuctuations 
of the entire system (e.g. in the context of biological 
systems this would be equivalent to large-scale protein 
breathing). The distributions of AErand foi" the eight 
clusters are shown in Fig. ^p. Strikingly, all points over- 
lapped into three regions that we colored in red, blue 
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FIG. 1. The quantum-structures network, (a) Nodes in the 
network represent structures with six excitable sites {N = 6) 
and E > 0.9. Links are placed if two structures are geo- 
metrically similar. The network layout is obtained via the 
Fruchterman-Reingold algorithm [17] which puts nodes with 
several neighbors in common close in space. To avoid over- 
crowding in the network layout, only 1/5 of the nodes have 
been represented in the picture. Node size is proportional to 
the number of links, (b) The probability distribution of effi- 
ciency loss upon random displacements for the eight clusters 
(see text for details). All eight distributions are shown where 
colors are chosen to highlight the presence of three classes 
(lines are a guide to the eye). Network nodes color coding fol- 
lows the class definitions: pair (red), inline (blue) and sparse 
(yellow) . 



and yellow. This coding split the network in homoge- 
neously colored parts, as shown in Fig. [l^. This was 
a-priori not obvious as the network could have shown a 
certain degree of color mixing, providing strong evidence 
that network topology and excitation dynamics are cor- 
related properties. Summarizing, the random displace- 
ments analysis provided evidence that (i) all structures 
within a cluster show similar response upon perturbation 
and (ii) three well-defined types of responses emerge from 
the eight clusters. 

These results suggested a characterization of the whole 
sample of efficient structures into three classes of similar 
quantum behavior. The classes are named "pair", "in- 
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FIG. 2. Structure-dynamics relationship for the (a) pair, (b) inhne and (c) sparse classes. On the left side of the figure 
structures belonging to the cluster containing the most efficient configuration of the class are overimposed (namely cluster 
number 1, 2 and 7 respectively). Characteristic structural motifs emerge from each of the different clusters (see Supp. Mat. 
for the superposition procedure). The exciton dynamics of the most efficient structure of each class is shown in the right part 
of the figure where the x and y-axis represent the time and the excitation probability, respectively. The relative frequencies to 
find structures with four, five or six active sites within each class are shown as insets (see text). Structural rendering was done 
with VMD 



line" and "sparse" because of their average geometrical 
properties. In fact, they respectively show couples of sites 
very close to each other, a compact arrangement around 
the input /output axis and a more sparse geometry. In 
terms of clusters, the pair class includes only one clus- 
ter (the most populated one, 73.3% of the whole sample) 
while the inline includes four clusters (17.6%) and the 
sparse only three (6.3%). 

The pair class is very robust against random displace- 
ments, with an average loss of efficiency of around 0.06 
(red data in Fig. [TJd). Interestingly, this class performed 
much better than the inline and sparse classes which 
showed an average loss of 0.14 and 0.10, respectively 
(with losses up to 0.3 for the former). These results in- 



dicated that the geometry in these latter classes is more 
correlated while in the pair class sites can be moved by 
small displacements in an independent fashion. However, 
robustness comes with a price: the pair class is generally 
slower. Taking as reference the time at which the maxi- 
mum in the output site is reached, the fastest processes 
for the pair, inline and sparse classes occur within 0.67, 
0.30 and 0.54 r, respectively. This pattern is also re- 
flected in the average speeds, with values of 0.92, 0.83 
and 0.86 r, respectively. 

In Fig. [2] several properties of the three classes are il- 
lustrated. On the left side of the picture structures be- 
longing to the cluster with highest efficiency within the 
class are overimposed on top of each other. This repre- 
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FIG. 3. Characterization of the pair class, (a) The most ef- 
ficient structure: pair and backbone sites are shown in dark 
and hght gray, respectively, (b) Efficiency loss upon pair re- 
moval as a function of the original efficiency E. Error bars are 
calculated according to the standard deviation, (c-d) Com- 
parison of the typical dynamics with and without the pairs. 
Coherent signal does not focus in the output site in the ab- 
sence of the pair sites, resulting in an efficiency loss AEpair 
of up to 0.32 (in the plot the dynamics of the most efficient 
structure is shown, with AEpair = 0.23). 



sentation allows a visual appreciation of the structural 
homogeneity within a cluster as well as of the diversity 
among clusters. In the right part of the picture, the exci- 
ton dynamics of the most efficient structure of the class 
is shown. In all three cases transport efficiency is larger 
than 0.97. However, the three dynamics differs substan- 
tially. In the pair class (right part of Fig. |2|i), two of 
the intermediate sites are successively excited with no 
active role of the remaining other sites (excluding input 
and output). Conversely, the other two classes show more 
complex patterns of excitation. These results provide ev- 
idence for a strong structure-dynamics relationship given 
that the final values of the efficiency are very similar in 
all three cases. Hence, structures differing from a geo- 
metrical point of view provide distinct exciton dynamics. 

As said, the pair class shows a prototypical modular 
structure. The first module is comprised of four sites 
including the input and output approximately lined up 
along the input/output axis, defining a backbone for the 
entire structure (light gray spheres in Fig. |3|i) . The sec- 
ond module is formed by the remaining two sites, the pair 
(dark gray spheres). Backbone sites are approximately 
equally spaced between input and output with typical 
inter-sites distances of around 0.60 — 0.64 tq. Pair sites 
instead are always very close to each other with a inter- 
site distance of around 0.25 ro (see Fig.[6]in Supp. Mat.). 



FIG. 4. The most efficient structures in the case N — A {E — 
0.922), AT = 6 (^ = 0.998) and AT = 8 (£; = 0.993). They 
are shown in light blue, gray and light green, respectively. 
The gray area highlights the presence of pairs for = 6 and 
AT = 8. 



The position of the pair within the box is more het- 
erogeneous than the backbone sites, i.e. their position in 
space changes between structures as indicated by the dis- 
perse cloud of sites in the structural overlaps of Fig. [2^. 
This property is connected to the exceptional robustness 
of this class upon random displacements of sites. The 
pair is usually found in a central position with a rela- 
tive displacement with respect to the input-output axis. 
Backbone and pair sites show a completely different dy- 
namical behavior. Systematic analysis of the distribution 
of the maximum excitations per site within the pair class 
(excluding input and output sites) revealed that back- 
bone and pair sites can be clearly divided into "active" 
and "inactive" exciton carriers, respectively (see Fig.[7|in 
Supp. Mat.). In fact, backbone sites present maximum 
excitation probabilities always larger than 0.25 while pair 
sites got never excited more than 0.075 in 99% of the 
cases. 

Strikingly, removal of the pair from the original struc- 
tures resulted in a systematic efficiency loss of up to 0.32, 
which is rather surprising given that pairs hardly serve 
as carrier of the excitation. Efficiency loss is particularly 
severe for the most efficient realizations due to the sensi- 
tivity of perfect constructive interference against pertur- 
bations (Fig. jsj)). Pair removal affects exciton transport 
in non-trivial ways as shown by the quantum dynamics of 
the most efficient structure before and after pair-removal 
(see Fig. [sj^-d). In the modified structure, we found un- 
balanced transport along backbone sites contrary to what 
was originally observed as well as an inability of the in- 
put site to promptly transmit the excitation to the closest 
backbone site. 

From a geometrical point of view the pair class provides 
attributes commonly observed in models with different 
number of sites. For N = A the ensemble of structures 
with E > 0.9 organized on a line with inter-site spacings 
very similar to the ones found in the backbone of the 
pair class. This is shown in Fig. [4] by superimposing the 
most efficient structures ioi N = A and A/" = 6 (light blue 
and gray spheres, respectively). However, the maximum 
efficiency ioi N = A is of only 0.922 while the presence 



5 



of pair sites lead to a maximum efficiency of 0.998 which 
is also the best value achieved within the whole sample 
of 100 millions structures. The general role of pairs is 
conffimed by the N = S case. In this model roughly half 
of the whole sample with E > 0.9 presented a modu- 
lar structure with four out of eight sites organized in 2 
pairs with the remaining sites perfectly overimposed to 
the backbone of = 6 (green spheres in Fig. [4|. Such 
geometry clearly represents a generalization of the pair 
class. Altogether, these results provided strong evidence 
that pair sites play an important structural role to tune 
up quantum efficiency, suggesting the idea that they rep- 
resent a general building strategy towards efficient trans- 
port in multi-body quantum systems. 

The presence of active and inactive carriers of the 
transport efficiency in the pair class pointed out a po- 
tentially intriguing scheme to rationalize the quantum 
dynamics of all classes. Taking 0.075 as a threshold to 
define a site as inactive, the pair class is characterized by 
a total of four active sites including the input and output 
(see probability distribution of active sites in the inset of 
Fig. |2|l). However, for the inline and sparse classes no 
clear separation into number of active sites was found 
(inset plots in panel b and c of Fig. [2j respectively). 
Moreover, when this concept was used for the analysis of 
the efficiency loss upon random displacements by group- 
ing structures according to the number of active sites, 
ambiguities raised between the five and six active sites 
groups due to strong overlaps (see Fig.jsjin Supp. Mat.). 
These results indicated that the active site concept alone 
is not sufficient to define structural groups with similar 
dynamics, reiterating the idea that advanced techniques 
for structural comparisons, like the one performed here, 
are necessary to unravel the connection between struc- 
ture and quantum dynamics. 

All the presented analysis was performed with per- 
fectly coherent dynamics, but it can be easily extended 
to models including dephasing (see Supp. Mat. for de- 
tails). We did so with different dephasing rates, ranging 
from 0.4/r to 2.0/r . Dephasing consistently decreased 
the efficiency with no specific distinction among the three 
classes, most importantly without affecting the response- 
to-noise classification identified with coherent dynamics 
(see Fig. [9] in Supp. Mat.). The purely destructive ef- 
fect of dephasing is consistent with the type of analysis 
performed, which focused on outstandingly fast and effi- 
cient transport made possible only by constructive inter- 
ference. 

Our work provided strong evidence for the emergence 
of non-trivial characteristic structural motifs leading to 
high quantum transport efficiencies. Specifically, the 
identification of site pairs not directly participating to 
the exciton dynamics results in a general strategy to both 
enhance quantum efficiency and to make structures ro- 
bust against geometric perturbations. Consequently, a 
design principle is presented, exploiting enhanced quan- 



tum transport in cases where perfect interferometric sta- 
bility is impossible. Gaining control on this problem is 
of paramount importance towards the rational design of 
technologies making use of constructive interference. 
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SUPPLEMENTARY MATERIAL 

Network Creation 

Network links are put according to a similarity param- 
eter. The similarity parameter S is calculated as: 



(3) 



where di is the distance between corresponding sites and 
i runs over all six sites. The sites are indistinguishable, 
thus all different permutations of the site labels need to 
be performed (in number of (A^ — 2)!, excluding input 
and output). In addition, symmetry along the in-out axis 
and an additional mirror symmetry has to be taken into 
consideration. The measure S between configurations A 
and B is thus calculated as follows: keeping fixed the set 
of labels for A, the labels of B are changed. For each of 
these (A^ — 2)! sets, 180 rotations of 2 degrees each of the 
configuration B are done around its in-out axis. For each 
rotation, a mirror refiection relative to the x — y = plane 
is also done. The final value of S is then the minimum 
value of Yl7=i among the (A^ — 2)! x 180 x 2 possible 
combinations of labels, rotation and mirror states. 

Only S'^ values below 0.0125 Tq are considered as links 
in the network. That means, the average di between sites 
of two "superimposed" configurations must be smaller 
than Y^O.0125 Tq ~ 0.11 tq. This value lies just above the 
tail of the pairwise distance distribution (Fig. [lO|, being 
only the most similar structures linked together. Lower 
values of the cut-off would generate a disconnected net- 
work, while values too close to the maximum of the dis- 
tribution would put links between structures that are not 
very similar. A double check with another cut-off value of 
S'^ < 0.0138 Tq was performed, giving results very similar 
to the ones shown in the main text (see Fig.|8|. Although 
the change from < 0.0125 to S'^ < 0.0138 seems 
negligible, it is worth noticing that it corresponds to an 
increase of the total number of links from 1.03 • 10^ to 
1.42 • 10^, i.e. a conspicuous 40%. 



Clusterization Procedure 

In order to obtain different homogeneous classes of 
configurations, the Markov cluster algorithm (MCL) was 
used [m |16] . This algorithm is based on the behavior of 
random walks on the network and consists of four steps: 



(i) start with the transition matrix A of the network, 
where each column is normalized to 1; (ii) compute A^; 
(Hi) take the pth power {p > 1) of every element of A^, 
normalizing afterward each column and (iv) go back to 
step ii. After some iterations MCL converges to Amcl, 
where only few entries for each column are non-zero (ex- 
actly only one non-zero entry per column) . These entries 
give the clusters. The parameter p is related to the granu- 
larity of the clustering process. High values of p generate 
several small homogeneous clusters. On the other hand, 
in the limit of p = 1, only one cluster is detected. In this 
work p = 1.4, giving a very good signal-to- noise ratio 
(see Fig. |5]Supp. Mat.). Smaller values of p give simi- 
lar results. For instance, p = 1.2 splits the network into 
two clusters: the one corresponding to the "pair" class 
described in the main text and the remaining two classes 
all together. This suggests that differences between the 
clusters in "inline" and "sparse" classes appear at a finer 
degree of granularity, while separation between these and 
the "pair" class is more evident, requiring a lower param- 
eter p to be resolved. 



Structural superposition 



Figures [2] and pT] are obtained as follows: for each clus- 
ter, the most connected structure is taken as reference 
and all the others are superimposed to that: for each of 
them the combination of labeling, rotation and mirror 
state which minimizes the similarity parameter S was 
taken (see Network Creation section) . In order to reduce 
noise, the coordinates of the sites were averaged with the 
ones from two other structures of the cluster taken at 
random. Structural rendering was done with VMD [18]. 



Dephasing Dynamics 

Loss of coherence in the system is introduced via the 
addition of an incoherent term to the dynamics of the 
quantum state. We consider an isotropic single body 
Markovian dephasing. The evolution of the quantum 
state in time is given by the Lindblad equation 



d 
di 



(4) 



where = ai^^ and ai^^ applies the Pauli matrix 
on the i-th site leaving the remaining unperturbed. 
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Cluster Ranking 

FIG. 5. Cluster populations. The populations of the clusters from 2 to 8 are fitted by a power law function (blue line). The 
first eight clusters represent cumulatively the 97% of the whole sample. The rest of the clusters, here fitted by a green line, are 
considered as noise. 




FIG. 6. (left) The most efficient structure of the pair class, (right) Inter-sites distances within the pair class and the N — 4 
efficient structures are depicted in gray and red, respectively. The distances dbb and ds show similar distributions in both cases 
(0.60 — 0.64 • ro for N = 6 and 0.59 — 0.61 • ro for A/" = 4), revealing a similar backbone geometry, while the intra pair distance 
dp is much smaller (0.25 • ro). 
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FIG. 7. Distribution of the maximum of the site excitation probability within time r (black, in and out sites are not included). 
The distribution of the pair class reveals two trends (red line): roughly half of the sites never gets excited more than 0.15 (the 
99% of which never more than 0.075), the other half has the maximum of excitation above 0.25. The sites belonging to the 
sparse (yellow) and inline classes (blue) stay in the middle of the distribution with no evident separations. 
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FIG. 8. Efficiency loss upon random displacements (see main text) according to the number of active sites (left) and the class 
partitioning proposed in this work (center and right), (left) Structures with 4,5 and 6 active sites are shown in red, yellow 
and blue respectively. This classification is not able to resolve the three different responses to noise presented in Fig. [TJ) of the 
main text (reproposed here and grouped by classes for comparison, center). Increasing the number of links of the network by 
a ^ 40% does not affect the behavior (right), suggesting our results are robust for different values of the similarity cutoff (see 
also Fig.IlO). 
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FIG. 9. Efficiency loss upon dephasing. (left) The three classes are shown in red, yellow and blue for the "pair", "sparse" and 
"inline" classes, respectively. Three different values of 7 were used (0.4, 1.3 and 2.0/t), here presented in solid, dashed and 
dotted lines, respectively. No evident difference in dephasing response subsisted between the classes, (right) Efficiency loss 
upon random displacement in the presence of dephasing (7 = 1.3/r). 




FIG. 10. Pairwise distance distribution between any- two structures for = 6 and E > 0.9. The chosen cutoff at 5* = 0.11 ro 
(5*^ = 0.0125 ro) is shown as a black vertical line. 
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FIG. 11. Superimposition of all structures belonging to clusters 1 to 4. The cluster number is colored according to the class 
of affiliation (red, blue and yellow for pair, inline and sparse classes, respectively). For clarity, structures are shown in three 
different orientations. 
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FIG. 11. (continued) 



Superimposition of all structures belonging to clusters 5 to 8. 



